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Abstract 

We study the effect of p° — 7 mixing in e + e~ — > it + tt~ and its relevance for the comparison of the square modulus 
of the pion from-factor \F^ e \s)\ 2 , as measured in e + e~ annihilation experiments, and |Fi T '(s)| 2 the corresponding 
quantity obtained after accounting for known isospin breaking effects by an isospin rotation from the r-decay spectra. 
After correcting the r data for the missing p — 7 mixing contribution, besides the other known isospin symmetry 
violating corrections, the irn 1=1 part of the hadronic vacuum polarization contribution to the muon g — 2 are fully 
compatible between r based and e + e~ based evaluations, r data thus confirm result obtained with e + e~ data. Our 
evaluation of the leading order vacuum polarization contribution, based on all e + e~ data including more recent BaBar 
and KLOE data, yields aJ; ad ' LO [e] = 690.75(4.72) x 1(T 10 (e+e" based), while including r data we find a); ad ' LO [e, r] = 
690.96(4.65) x 10 -10 (e+e~+r based). This backs the ~3 a deviation between a c * pcrimcnt and a* eory . For the r di-pion 
branching fraction we find = 25.20 ± 0.0.17 ± 0.28 from e+e~+CVC, while B^o = 25.34 ± 0.0.06 ± 0.08 is 

evaluated directly from the r spectra. 

Key words: 7 — p mixing, p-meson properties, e + e _ -annihilation, r-decay pion form factor, muon anomalous magnetic 
moment. 

PACS: 13. 66. Be, 13.35. Dx 14.60.Ef 



1. Introduction 



Isovector data for the pion form factor obtained from hadronic r-decay spectra can be compared with the 
mixed isovector-isoscalar data measured in the e + e~ channel by means of theory input [1]. In particular, we 
need some model in order to be able to disentangle p — u mixing as well as other isospin breaking (IB) effects. 
The general problem in confronting measured quantities like \F^ T \s)[I — 1]| 2 and |.Fi e ^(s)| 2 = \F^ e \s)[I = 
1] + Fi e \s)[I = 0]| 2 is the fact that the latter object is subject to quantum interference between the two 
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amplitudes and in general may not be well approximated by \Fi e \s)[I = 1]| 2 + \F^ e \s)[I — 0]| 2 . Without a 
specific model for the complex amplitudes one cannot get the precise relationship. 

Commonly, pion form factors measured in the neutral channel in e + e~ — > tt + it~ and in the charged 
channel in r~ — > v t tt~tt decay (or its charge conjugate) are parametrized by an extended Gounaris-Sakurai 
(GS) formula 

_ BW p G ( s 770) ( S ) • (l + ^BW U («)) + P BW« ( S 1450) ( S ) + 7 BWf 17m) (s) 
F7r{s) 1 + + 7 ' (1) 

which results as a sum of mixing isovector states, each described by a Breit-Wigner (BW) type of amplitude. 
The pion form factor is related to the corresponding cross section by 1 

o(e+e- 7T+7T-) = ^ ^ |Fi e )( S )| 2 = ^ v (s), (2) 

6 S S 



for point-like pions F n (s) = 1, where /3 n is the pion velocity in the cm. frame: /3 n = yl — 4m 2 /s. The 
spectral function Uj(s) is related to the form factor by 

«i(*) = ^l4°WI 2 ; (i = 0,-)+>(e,r) , (3) 

for the neutral (0) e+e~- and charged (-) r-channel. The spectral function V-(s) can be measured very 
precisely in T-decay: 

ldr _ 6\V ud \ 2 S EW B e ( s \ ( 2s\ 

with m T = (1776.84 ± 0.17) McV the r mass, \V ud \ = 0.9418 ± 0.00019 the CKM matrix clement, B e = 
(17.818 ± 0.032) % the electron branching fraction, B^^ = (25.51 ± 0.09) % the di-pion branching fraction 
and Sew = 1-0235 ± 0.0003 the short distance electroweak correction. 
Note that a single standard Breit-Wigner resonance yields 

2 _36 T(p^e+e-) s sT 2 p 
1 w{ > l o? PlT{p -»■ TT+TT-) Mj (s - M 2 ) 2 + M 2 r 2 ' 

Denoting the 7 — p transition coupling by eM 2 / g p the branching fraction at resonance reads 

r(p^e+e-) a 2 /^\ 2 f M.V 
np r(p^7r+7r-) 36 J 

with (3 p = (3 v (s — M 2 ). In the case of complete p dominance g p — g p7rw 2 . 

The GS formula (1) also describes the charged isovector channel provided 6 = 0, since there is no charged 
version of the to. In the neutral channel the GS formula does not fully include p° — 7 mixing , which is known 
since the early 1960's, when the p had been discovered. A direct consequence of p° — 7 mixing is the vector 
meson dominance (VMD) model characterized by an effective Lagrangian [4] 

eM 2 

C 1P = P -P»A* . (5) 

9p 



1 QED corrections to e+e~ — > w+n~ have been summarized in [2,3]. They will not be considered in the following and we 
assume them to be taken into account in the extraction of the form factor from the experiments. 

2 At resonance the single BW pion form factor is given by 

and for PDG values of the parameters yields | F„ (M 2 ) | 2 ra 39 a reasonable value (see below) . 
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However, this form does not preserve electromagnetic gauge invariance and the photon would acquire a mass 
unless we add a photon mass counterterm to the Lagrangian which is fine tuned appropriately. The pion 
form factor here takes the form 

and the condition of electromagnetic current conservation ^(0) = 1 is satisfied only if g p7r7r — g pi which 
is called universality condition and corresponds to complete p dominance. In fact electromagnetic gauge 
invariance can be implemented by writing the effective VMD Lagrangian in the form [5] 

z Up 

in terms of the field strength tensors. As it satisfies gauge invariance, the form factor calculated here reads 

and satisfies the current conservation condition F ff (0) = 1 in any case, irrespective of the universality 
constraint g p7nT = g p (for a recent discussion also see [6]). Obviously, this simple model is not able to 
describe the pion form factor measured in e+e~ — > tt + it~ at low energies, unless we take into account 
energy dependent finite widths effects of the p as it is done in the GS model [7] 3 . The energy dependence 
of the p-width has to reflect the off-shell p* — > ttt: process. So we have to model effectively a "rho-pion- 
photon" system, discarding the co and its mixing with the p, which is well understood and will be taken 
into consideration in a second step. Our focus here is to work out the difference in the relation between the 
charged channel and the neutral channel, which results from the p° — 7 mixing . The latter like p° —to mixing, 
has no counterpart in the charged channel. The purpose of this study is to understand better the discrepancy 
between r and e+e~ di-pion spectra, which has been clearly established in [16] under the assumption that all 
possible IB corrections were accounted for. More recent data form Belle [17] and KLOE [18,19], and applying 
improved IB corrections, confirmed a significant discrepancy [20]. Although the new tt + tt~ spectrum from 
BaBar [21], measured via the radiative return mechanism, is closer to the corresponding spectra obtained 
from r-decays, a discrepancy persists. 



2. A p — 7 mixing model and related self-energy effects 



As already said, the VMD ansatz has to be replaced by a more realistic model which must take into 
account 

- the finite p-width, related to its decay p — > , 

- the p — 7 mixing, which leads to non-diagonal propagation of the p — 7 system, and 

- the p — uj mixing, which we will consider in a second step. 

This has to be implemented in an appropriate effective field theory (EFT). In a first step we consider the 
interaction of the p with the pions together with their electromagnetic interaction, assuming the pions to 
by point- like (scalar QED). As suggested long ago by Sakurai [22], the p may be treated as a massive gauge 
boson. The effective Lagrangian thus reads 

C = C 1P + ; £ v = D p ir + D +t *ir~ - m*-K + 'K~ ; = d p - ie A p - ig p7T7T p p . (9) 
The corresponding Feynman rules in momentum space are 

3 Other models have been reviewed and investigated recently with emphasis on p — ui mixing in Ref. [8]. Frequently used 
descriptions of the low energy irir form factor include the ChPT-based Guerrero-Pich formulation [9], the Leutwyler-Colangclo 
approach [10], the resonance Lagrangian approach [11] (see e.g. [12]) or the related Hidden Local Symmetry (HLS) model as 
applied in [13,14], and the phenomenological Kuhn-Santamaria (KS) model [15]. As we will see our model is closely related to 
the GS model, and we adopt the latter for comparisons and fits. 
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A^rnr = -ieip + p'Y ; p^mr = -ig^n (p + p'Y 

A fJ, A l, TTiT = 2ie 2 g^ ; ffp u mr = 2\g 2 p ^g^ 

A" (Twit = 2ieg p ^g^ ; A»p v = -i e/g p (p 2 g" v - p»p v ) . 
The model should be understood as a simplified version of the better justified effective resonance Lagrangian 
approach [11], which extends the chiral structure of low energy QCD (chiral perturbation theory) to include 
spin 1 resonances in a consistent way. A variant is the HLS model, which in the same context has been 
applied to investigate the (p,ui,(f>) mixing effects in [13]. Actually in [13] too, V — 7 (V — p, lj, cf>) mixing 
amplitudes have been included (more on that below). The main difference to the GS model is that we take 
our EFT Lagrangian serious in the sense that we include all relevant contributions to e + e~ — > tt + t:~ , while 
in the GS model some of the contributions have been neglected. In fact the GS model is incomplete in the 
sense of a quantum field theory. 

In sQED the contribution of a pion loop to the photon vacuum polarization is given diagrammatically by 



-in^W(g) 



and one then obtains the bare 7 — p transverse self-energy functions 



n 77 = 7^/(9 2 ) . n 7P = ^/(<7 2 ) and U pp = %\ f( q 2 ) , (10) 



where 

f(q 2 ) = q 2 h(q 2 )= ^B (m^m„;q 2 )(q 2 -4ml)-4A (m„)-4ml + ^ , (11) 

in terms of the standard scalar one-loop integrals Ao(m) and Bq^tti, m; s ) [23]. Explicitly 4 , in the MS scheme 
(p the MS renormalization scale) 

h(q 2 ) = f(q 2 )/q 2 = 2/3 + 2 (1 - y) - 2 (1 - y) 2 G(y) + In i£ , (12) 

where y — Am 2 / s and G{y) — ^g- (In pffi — i 7r), for q 2 > 4 m 2 . Note that all components of the (7, p) 2x2 
matrix propagator are proportional to the same function f(q 2 )- The renormalization conditions are such 
that the matrix is diagonal and of residue unity at the photon pole q 2 = and at the p resonance s = M 2 , 
hence the renormalized self-energies read (see e.g. [24]) 



n? 7 V) 


= n 77 (? 2 )- 


q 2 n' 11 (0)=q 2 ri^(q 2 ) 


(13) 


n^pV) 


= n 7P ( 9 2 )- 


^Ren 7P (A/ p 2 ) 


(14) 


n p c p n (<z 2 ) 


= n pp ( 9 2 )- 


ReII pp (A/ 2 ) - (q 2 - M 2 ] 


,Re^(M 2 ) (15) 


where IT 77 (0) 


= n 7P (o) = 


n pp (o) = and n; 7 ( g 2 ) 


= H 17 (q 2 )/q 2 , has been used. Note, that the tree level 



mixing term in the Lagrangian contributes to the bare 7/3 self-energy as LI 7P ^ = q 2 (e/g p ), which does not 



4 The standard Gounaris-Sakurai parametrization differs from our sQED model and utilizes h(q 2 ) = — 2 (1 — y) 2 G(y) which 
for q 2 — ¥ behaves as h(q 2 ) — > 2y = fira^/q 2 i.e. q 2 h(q 2 ) — > Sm 2 , which in D 77 represents a non- vanishing photon mass. While 
the constant terms in (12) drop out by renormalization, the +2(1 — y) term is required by electromagnetic gauge invariance 
and in fact renders h(q 2 ) — ¥ const, regular in the static limit. 
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affect the renormalized self energies, however. In particular, <5U 



rcn 
IP 



function renormalization reads Z p = 1/(1+ -tjjM 5 = M 2 )) with 



an 



As 



M 2 P ) 



9p7T7T 

48 7T 2 



8/3 - p 2 p 



1 + (3 -^JK ln 



1 + & 



q 2 — ~ -wi M l— = 0- The p wave 



(16) 



Numerically, Z p ~ 1.1289 at p, = in,. 

It is crucial to observe that vacuum polarization effects affect mass renormalization of the p as well as 7 - 
p mixing, in spite of the fact that photon vacuum polarization has to be subtracted in the definition of F^. 
In other words, in sQED we would still have F 7T (s) = 1 in (2) while vacuum polarization is absorbed into a 
running fine structure constant 



a — > a(s) 



1 + n;™ n (s) 



(17) 



which mean that in calculating F 7T (s) we have to multiply the result by 1 + H™ n (s). 
A convenient representation of U p c p n is given by 



7TM p ^ p 



s (h(s) - Re h{M 2 p )) - (s - M 2 p ) M 2 p Re 



dh 
As~ 



--Ml 



with 



s^(s) = 3y-l-3y(l-y)G(y) . 



In particular 5 : 



n ren (o) = MpTp M 2 — 

pp( ' tt/33 m P As 



Mp£p 



(3y p -l-3y p /3 2 p G(y p )) ; y p = 4m 2 JM 2 



(18) 



(19) 



(20) 



Without mixing, pion production mediated by the p resonance, yields the GS type pion form factor, nor- 
malized to ^(0) = 1, 



-M 2 +U^(0) 



s-M 2 + U p f( S ) ■ 
The renormalized mixing self-energy may be written in a form 



( ' ~ 48tT 2 



{s (h(s)-Reh(M 2 ))} . 



(21) 



(22) 



Note that while the inverse propagator matrix is diagonal at the two propagator poles, off the poles it is 
not diagonal. This is the main effect we are going to discuss now 6 . 

The propagators are obtained by inverting the symmetric 2x2 self energy matrix 



5 In contrast to sQED, in the standard GS formula 



dh 



2(1 -y) 2 G(y) ; s — (s) = y - 1 - 3y(l - y) G(y) , 
ds 



which is singular for s — > 0. The first term in (18) in this case yields a finite contribution sh(s) — > 8m^ and thus 



M T 



8ml 



p As 



s=Ml 



(3y p -l-3y p fiG(y p )) 



i.e., ripp n (0) = — dT p M p is actually not modified, in spite of lacking manifest gauge invariance. 

6 These effects are very similar to Z° — 7 mixing [25] which has been investigated theoretically as well as experimentally at 
LEP with high precision (sec Rcfs. [26,27] and references therein). Typically, these effects at low energy or near the Z pole are 



£_i = | <7 2 + n 77 ( g 2 ) n 7P ( g 2 ) , 3 



n 7P (^) ^-M; + n pp (^) 



with the result: 



Z? 77 — 



«2 + n ^ n ^ (,?2) <7 2 + n 77 (g 2 ) 



-n 7P (g 2 ) ^ -n 7P (g 2 ) 



w (g 2 + II 77 (g 2 ))(g 2 - M p 2 + II pp (g 2 )) _ n 2 p (g 2 ) " (g 2 + II 77 (g 2 ))(g 2 _ M 2 + n w ( ? 2)) 

^ = f-Mt + n^)-^^ * ^MJTIW) • (24) 

These expressions sum correctly all the irreducible self-energy bubbles 7 . The approximations indicated are 
the one-loop results. The extra terms are higher order contributions and are particularly relevant near the 
resonance, characterized by the location sp of the pole of the propagator, which is given by the zero of the 
inverse propagator: 



n 2 (sp) 

s P - m 2 „ - n p o p o (s P ) 7P „ . : = , (25) 

H sp -1I 77 (Sp) 

with sp = M p0 complex. The usual (no mixing) considerations in determining the physical mass and width 
of a resonance remain true if we denote the self-energies by II y (V = p°,p ± ) with 

n p± (/,-•■) = n p+p - (p 2 ,...) 

and 

Thus the location of the pole may be written as 

M 2 -m 2 +n(M 2 ,m 2 ,---) = 0, (26) 
for both the p ± and the p°, where 

M 2 p = (q 2 ) pole = M 2 p -iM p T p 

is characterized by mass and width of the p 8 . Note that the imaginary part of the self-energy function is 
energy dependent, which implies an energy dependent width, of course with the correct phase-space behavior 
of p — > 7T7T decay. 

How do off-diagonal elements of the 7 — p propagator affect the line-shape of the pi We assume we know 
the p mass M p and the p width T p for the unmixed p as it is seen e.g. in the isovector r decay spectra, i.e. 
we compare the result with a charged p^- assuming equal mass and width. We therefore compare result first 
with the Belle data [17]. Of course our model does not fit the data, because a more sophisticated extended 



expected to be small because of the smallness of the electromagnetic fine structure constant, and at the Z resonance, because 
of the large mass and very small width of the Z° boson. In case of the p witch's mass lies not very far above the hadronic 
7T7r-threshold (which is very low by the fact that pions are quasi Nambu-Goldstonc bosons) and due to the relatively large 
(hadronic) width we expect corresponding mixing effects to be much more relevant. In the charged channel, in principle, there 
is W ± - p ± mixing, with some effective W + p + h.c. coupling term. However, this produces a negligible effect because the 
W propagator pole is far away from the p propagator pole and from the energy range of interest. In fact the mixing matrix, 
diagonalized at the p-pole, remains essentially diagonal in the whole range of interest (< 2 GeV). 

7 It is of curse well known that this Dyson summation is crucial for a proper description of the particle/resonance structure in 
particular near the poles, where naive perturbation theory in any case breaks down. 

8 For 7 - Z° mixing in the electroweak Standard Model explicit results up to two loops have been worked out in [28] . 



G 



40.0 



30.0 



20.0 



10.0- 




ielle 08 
no p" (y=0) 
p only (y=/?=0) 



400.0 800.0 1200.0 

E (MeV) 



1600.0 



Fig. 1. GS fits of the Belle data and the effects of including higher states p' and p" at fixed M p and T p . Doubts are in order 
whether the higher resonances really affect the p resonance in the way suggested by the commonly adopted GS paramctrization. 

Gounaris-Sakurai model (1) has been used to extract the p parameters. If we switch off the contributions 
from p' and p" by setting 7 = and (3 — 7 = we observe a substantial change in F n (s) as illustrated in 
Fig. 1. Note that in an EFT one would expect the heavier states to decouple, while in the GS type modeling 
the low energy tail is normalized away by the 1 + j3 + 7 normalization factor. In field theory in place of this 
normalization a factor s/M r 2 es would imply automatic decoupling. But that is not the way mass and width 
of the p are determined usually. Evidently, in the GS model, in the p-region the higher resonances serve as a 
continuum background without which good fits in general are not possible. So if we stick with our simplified 
model we cannot expect to get a good representation of the data without corresponding extensions. On the 
other hand, the simplified model allows us to work out more clearly the effect of 7 - p mixing. Above we 
have diagonalized the mixing propagator matrix at the poles, this allows to make precise the meaning of 
mass and width of the heavy unstable state. This is achieved if renormalized mixing self-energy is given by 



n 7P ren ( q 2 ) = n 7P ( g 2 ) - n 7P (o) - (i? e n 7P (M 2 ) - n 7P (o)) . 

This can be achieved by two subsequent transformations of the bare fields: 
i) Infinitesimal (perturbative) rotation 



(27) 




diagonalizing the mass matrix at one-loop (n+l-loop) order given that the mass matrix has been diago- 
nalized at tree (n-loop) level, 
ii) Upper diagonal matrix wave function renormalization inducing a kinetic mixing term (this cannot be done 
by an orthogonal transformation) 




which allows to normalize the residues to one for the 7- and p-propagator, respectively, and to shift to 
zero the mixing propagator at the p-pole. 
Thus the relationship between the bare and the renormalized (LSZ) fields is (expanded to linear order) 



A b = ^Z, ( A r - (A p + A ) Pr 
Pb = \fZp~Pr + A A r , 



(28) 
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generalizing the usual multiplicative field renormalization represented by the first term for both fields. The 
counter-terms Aq and A p are determined by the condition (27) 



A : 



n 7P (Q) 
mi 

ReXl lp {M 2 p ) - n 7P (0) 
M 2 p 



(29) 



For our model Ao = and A p = e/g p to leading order. The field transformations of course induce mixing 
counter terms at the vertices, which are absorbed into the definition of the physical couplings. In principle, 
this non-symmetric transformation only affects the bookkeeping such that the propagator pole structure 
becomes obvious. It does not change the value of the functional integral i.e. the mixing counter terms cancel 
in the interior of Feynman diagrams, unless the photon and/or the rho are involved as external fields (states). 

As a consequence of the diagonalization the physical p acquires a direct coupling to the electron: starting 
as usual from the bare Lagrangian 

£qed = ^e7 M (^ - i e b A bp ) V> e (30) 
we obtain 

£qed = ^eT^M - ie^4 M + igpeep^ipe (31) 

with g pee — e (A p + Ao), where in our case Ao = 0. 

The e + e~ — > tt + tt~ matrix element in sQED is given by 

M = -ie 2 v^u{ Pl ~p 2 ) p F n {q 2 ) (32) 

with F 7r (q 2 ) — 1. In our extended VMD model we have the four terms shown in Fig. 2 and thus 




Fig. 2. Diagrams contributing to the process e + e 



i^Ti-(s) OC 6 Dyy -)- Cf7 P 7r7r D '-y p g pee^-F) p^ Q peeQ p-KTT ^ pp ) 

where the first term properly normalized must be unity 9 . Thus 

^7r(s) = Dryy ~\~ 6 (^/pTTTT Q P^e) Q peeQ pTZTT D ppj / Dryy^ 



(33) 



Note the sign of the induced coupling g pee in (31), which leads to the signs as given in (33). Typical couplings 
read gp^barc = 5.8935, g p ^ IC n = 6.1559, g pee = 0.018149 and x = g^/Op = 1.15128. 

Real parts and moduli of the individual terms normalized to the sQED photon exchange term are displayed 
in Fig. 3. 

An improved theory of the pion form factor has been developed in [10]. One of the key ingredients in 
this approach is the strong interaction phase shift Sl(s) of tttt (re)scattering in the final state. In Fig. 4 
we compare the phase of F 7r (s) in our model with the one obtained by solving the Roy equation with tttt- 
scattering data as input. We notice that the agreement is surprisingly good up to about 1 GeV. It is not 
difficult to replace our phase by the more precise exact one. 

We note that the precise s-dependence of the effective p-width is obtained by evaluating the imaginary 
part of the p self-energy: 



9 Note that the conserved vector current (CVC) condition -F,r(0) = 1 in our model is given and saturated by the sQED term 
Dj^y alone, while in the GS model -FV(O) = 1 is imposed by force on the term D pp , the only one present in the GS case. 
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0.0 200.0 400.0 600.0 800.0 1000.0 
E (MeV) 

Fig. 4. The phase of F 1T (E) as a function of the c.m. energy E. We compare the result of the elaborate Roy equation analysis 
of Ref. [10] with the one due to the sQED pion-loop. The solution of the Roy equation depends on the normalization at a high 
energy point (typically 1 GeV). In our calculation we could adjust it by varying the coupling g pwlr . 
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!mlf ; .... = ^^ s = M,r.4.M - C-!4) 

48 7T 



which yields 



2 2 



P 

In our model, in the given approximation, the on p-mass-shell form factor reads 



F T (M 2 ) = l-i^f^^ , (36) 

6 1 p 



and the square modulus may be written as 

36 T ee 

with 



F.(M 2 )| 2 = l + §-^ , (37) 



Tpee = \ 9j ^ M P 0r 9pee = ^l2vT pee /M p . (38) 

It is interesting to note that the GS formula (21) does not involve T pee in any direct way, since the nor- 
malization is fixed by applying an overall factor 1 + dT p /M p = 1 — n pp n (0)/M 2 to enforce i^(0) = 1. The 
leptonic width is then given by 



pee g p 



r£ = n ; e (i + dv p /M p y . (39) 



p 

In the CMD-2 fit 1 + dT p /M p ~ 1.089. 

The result for ^(s)! 2 (using mass and width as before) is displayed in Fig. 5. We compare the results 
obtained when p - 7 mixing is properly taken into account with the one obtained by ignoring mixing and 
with a GS fit with just the p taken into account. At first look, the results agree fairly well but do not fit the 
Belle data as expected if we do not include the higher resonances. 

A detailed comparison, in terms of the ratio 

r (a)- |F - (S)|2 (40) 

shown in Fig. 6, however reveals substantial differences and proves the relevance of the mixing. We also 
plotted the same ratio for the; 1=1 part of the GS fit, which exhibits a similar behavior as the true: F v (s). 
This is not really surprising as one fits the same data just in a different way, i.e. with different parameters 
for the p. This mixing affects, however, the relationship to the r channel, which does not exhibit this effect. 
Of course at higher energies, not to far above the p, it is not known whether the simple EFT model can 
be trusted. Note that dropping the II 2 p terms (approximation indicated in (24)) in the Dyson resummed 
propagators does not affect the result. We have checked that lo — 7, <p — 7 or T(45) — 7 mixing effects are 
tiny away from the resonances and thus should not affect the interpretation of radiative return spectra as 
measured at KLOE and BaBar. 

We have to compare the model with the 1 = 1 part of the e + e~-data. To this end we may take the CMD-2 
fit of the CMD-2 data [29] and set the mixing parameter S = as illustrated in Fig. 7. In this way we obtain 
the isovector part of the square of the pion form factor \F!f\l = l](s)| 2 . 

In order to compare Fi T \s) extracted from r-decay spectra with Fi e \s) measured in e + e~-annihilation 
we have to apply isospin breaking corrections as investigated in [30] and [31] (see also [16,32,20]): 
- Mass shift: the Cottingham formula, which allows us for a rather precise calculation of the electromagnetic 

pion mass shift Sm^ = m^± — m^o ~ 4.6 MeV, suggest the relation Am 2 ~ AM 2 , which yields a shift 

1 A 2 

between changed an neutral p by 6M p = M p ± - M p o ~ \ -j^f = 0.814 MeV. 
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Fig. 5. Effect of 7 - p mixing based on the simple EFT model (9). Parameters: M p = 775.5 MeV, T p = 143.85 MeV, 
B[(p -> ee)/(p -> tttt)] = 4.67 X 10~ 5 , e = 0.302822, g p7V7v = 5.92, g pee = 0.01826. The crucial point is the difference between 
"p no mixing" [1st plus 3rd term of (33)] and "p — 7 mixing included" [all terms of (33)]. The interference with the mixing 
term lowers the form factor above the y/s ~ M p , an effect not present in the charged (pure 7 = 1) channel. For comparison the 
GS fit with switched off p' and p" and the "p only" [3rd term of (33) only] are shown. 



Width shift: The kinematic shift form pion and rho mass differences is ST p = T p ± — T p a = (/3 3 ± M p ± 
P^o M p o) = 1.300 MeV. 

(m 7r ±)-m 2 



Pion velocities are fi^ = \ 1 — 



(2m +) 2 



and fi w ± = V 1 - 



(m i+m^o) 2 



for neutral and 



charged channel, respectively. Their on-resonance values read j3 p ± = /3„ t ±(s = M 2 ± ) and /3 p o = f3„fl(s 



M 2 pQ ) 



In the charged channel (r-dccay) the appropriate phase-space for the n tv system, replacing the tt + it 
one, has to be considered. For the energy dependent width one has T p ± (s) = T p ± 



is charge indepen- 



We have made use of the fact that the strong coupling factor -^ 2L = p , yi — - 

b f b> 487r MpQ M ^ ± fj.i ± 

dent. Note that all of these corrections represent corrections in Fq(s)/F-(s) the ratio between neutral and 
charged channel F^'s. 

- Electromagnetic corrections Gem(s) as calculated in [30,31]. Specifically, we will apply the correction 
given by [30], since the ones given in [31] differ quite a lot for reasons we have not yet understood. It does 
not affect the main conclusion of our analysis. 

In total a correction 10 



v (s) = Rm(s) v-(s) ; Rib(s) 



1 



Gem(s) p±{s) 



F-(s) 



(41) 



10 If we would not include the p — 7 mixing in Fq(s) the correction formula would read 

«o(s) = r P f(s) Ri B (s) v-(s) . 
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Fig. 6. a) Ratio of the full \F^(E)\ 2 in units of the same quantity omitting the mixing term together with a standard GS fit 
with PDG parameters, b) The same mechanism scaled up by the branching fraction Yy /T{V —> tttt) for V = u> and <j>. In the 
7T7r channel the effects for resonaces V ^ p are tiny if not very close to resonance. 
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Fig. 7. CMD-2 data for |.FV| 2 in p — u) region together with Gounaris-Sakurai fit. Left before subtraction right after subtraction 
of the u>. 



has to be applied in the relation between the spectral functions. Final state radiation correction FSR(s) and 
vacuum polarization effects we have been subtracted from all e + e~-data. 

In Fig. 8 we illustrate the consequence of p — 7 mixing. After applying the correction (for our set of 
parameters, which is not far from standard GS fit parameters) the consistency of r and e + e _ data is 
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Fig. 8. \F n (E)\ 2 ratio r vs. e+e - 1=1 (CMD-2 GS fit): a) uncorrected for p — 7 mixing, and b) the same after correcting for 
it. The correction factor is given by the solid curve of Fig. 6 a). Lower panel: e+e — energy scan data [left] and e+e — radiative 
return data [right]. The GS fit chosen as a reference is represented by the full line in the right Fig. 7. The 1=1 part for the 
e + e — sets is obtained by subtracting the difference of the two curves shown in Fig. 7. The choice of the particular reference is 
of course ambiguous. 

dramatically improved. However, substantial differences of different measurements remain as a problem. 

How does the new correction affect the evaluation of the liadronic contribution to the anomalous magnetic 
moment of the muon? To lowest order in terms of e + e _ -data, represented by R(s), we have 

00 

«T' LO M = ^ / dsi^OO^ , (42) 

with the well-known kernel K(s) and 

RW(s) = (3 S a™)/W( S )) = 'Sv ( S ) . 

Note that the p — 7 interference is included in the measured e + e~-data, and so is its contribution to ajj ad . In 
fact a^ ad is intrinsic an e + e~ -based "observable" (neutral current channel). If we want to use the r spectral 
function 11 isospin breaking corrections to the r data must be applied: traditionally v-(s) —¥ vq(s) — 
Rib(s) v -( s )- For our comparison we switch off the 1=0 (p — u) mixing part) part from the e + e~ data, which 
we do not include in Rib- In addition we have to account for the missing p — 7 mixing in the r spectra. 
The results for the 1=1 part of ajj ad is given in the Tab. 1. Remarkably, the contributions to ajj ad from the 
e + e _ -data on the one hand and from the r data on the other hand agree surprisingly well after including 
the mixing. The p — 7 mixing correction has been applied for the p part only (p 1 and p" subtracted and 
added back with the help of the GS fit from Belle). We do not attempt here a complete analysis of all r 



11 In the exact SU(2) (isospin) limit we would have the CVC relation ti_ (s) = vq(s) and the t data would give the same 
contribution as the e+e _ data. 



13 



Table 1 

Isovector (1=1) contribution to a p ad X 10 10 from the range [0.592 - 0.975] GeV from selected experiments. First entry: results 
from T-data after standard isospin breaking (IB) corrections. Second entry: results from r-data after applying in addition the 
p-7 mixing corrections r p7 (s), with fitted values for M P ,T P and F pee [M p = 775.65 MeV, T p = 149.99 MeV,B[(p -» ee)/(p -> 
7T7r)] = 4.10 X 10~°]. Third entry: as second one, however using PDG values for the p parameters [M p = 775.5 McV,r p = 
143.85 MeV,B[(p — > ee)/(p — > nn)] = 4.67 X 10~ 5 ] (for illustration of sensitivity to the choice of parameters). For the p — uj 
mixing we subtracted 2.67 X 10~ 10 . Errors are statistical, systematic, isospin breaking and p — 7 mixing, assuming a 10% 
uncertainty for the latter. Final state radiation is not included. 



Data 


standard IB corrections 


incl. p — 7 mixing 


same using GS fit parameters for the p 


ALEPH 1997 [33] 
ALEPH 2005 [34] 
OPAL 1999 [35] 
CLEO 2000 [36] 
BELLE 2008 [17] 


390.75(2.69)(1.97)(1.45) 
388.74(4.05)(2.10)(1.45) 
380.25(7.36)(5.13)(1.45) 
391.59(4.16)(6.81)(1.45) 
394.67(0. 53)(3.66)(1. 45) 


385.63(2.65)(1.94)(1.43)(0.50) 
383.54(4.00)(2.07)(1.43)(0.50) 
375.39(7.27)(5.06)(1.43)(0.50) 
386.61(4.11)(6.72)(1.43)(0.50) 
389.62(0.53)(3.66)(1.43)(0.50) 


385.05(4.01)(2.08)(1.44)(0.36) 

388.10(4.12)(6.75)(1.44)(0.36) 
391.06(0.52)(3.63)(1.44)(0.36) 


average 


391.06(1.42)(1.47)(1.45) 


385.96(1.40)(1.45)(1.43)(0.50) 




CMD-2 2006 [29] 
SND 2006 [37] 
KLOE 2008 [18] 
KLOE 2010 [19] 
BABAR 2009 [21] 




386.34(2. 26)(2.65) 
383.99(1.40)(4.99) 
380.24(0. 34)(3.27) 
377.35(0. 71)(3.50) 
389.35(0. 37)(2.00) 




average 




385.12(0. 87)(2. 18) 




all e+e - data 
e+e~~ + t 




385.21(0. 18)(1. 54) 
385.42 (0.53)(1.21) 





data. The p — 7 mixing contribution, which implies a moderate positive interference below the p resonance 
and a stronger negative one above the p resonance, shift the r data to lie perfectly within the ballpark of 
the e + e~ data 12 . 

Another important quantity which can be directly measured is the branching fraction = Y(t — > 

v t ttit )/T t . This "r-observable" can be evaluated in terms of the 1=1 part of the e+e~ — > ir + n~ cross 
section, after taking into account the IB correction v (s) — > V-(s) = v (s)/ Rib(s), 



_B cy c _ 



25 , E w-B e |K(i| 



mt 



I 



dsR 



(0) 



■ W 1 



mt 



1 + 



2s 



mi 



Rm(s) 



(43) 



where here we also have to "undo" the p — 7 mixing which is absent in the charged isovector channel. 
Results are given in Tab. 2. The shift by the p — 7 mixing is +0.62 % to which we assign an error of 10 % 13 . 
For the direct r branching fractions the first error is statistical the second systematic. For e + e~+CVC the 
first error is experimental the second error includes uncertainties of the IB correction +0.06 from the new 
mixing effect. Remaining problems seem to be experimental, there are significant differences in the spectral 
functions from different experiments (see Fig. 8). 



12 Note that the moderate shifts of rho masses and widths [32] which diminish somewhat the discrepancy are included in the 
IB corrections, as detailed above. 

13 Averages given in [20,38] are 25.42 ± 0.10 % for r and 24.78 ± 0.28 % for e+e~ + CVC. Adding the new correction we get 
25.40 ± 0.28 ± 0.06 % for the latter, in perfect agreement with the r result. The BaBar di-pion spectral function, not included 
in the previous numbers, yields 25.15 ± 0.18 ± 0.22 % or 25.77 ± 0.18 ± 0.28 including the new correction. Results differ 
slightly from ours because we apply slightly different IB corrections. 
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Table 2 

Calculated branching fractions in % from selected experiments. Experimental data completed down to threshold and up to m T 
by corresponding world averages where necessary. CMD-2 is a combination of 2003 and 2006 data. The experimental world 
average of direct branching fractions is B CV C = 25.51 ± 0.09 % [39]. 



t experiments 


B^ [%] 


e + e experiments 




ALEPH 1997 


25.27 ± 0.17 ± 0.13 


CMD-2 2006 


25.40 ± 0.21 ± 0.28 


ALEPH 2005 


25.40 ± 0.10 ± 0.09 


SND 2006 


25.09 ± 0.30 ± 0.28 


OPAL 1999 


25.17 ± 0.17 ± 0.29 


KLOE 2008 


24.82 ± 0.29 ± 0.28 


CLEO 2000 


25.28 ± 0.12 ± 0.42 


KLOE 2010 


24.65 ± 0.29 ± 0.28 


Belle 2008 


25.40 ± 0.01 ± 0.39 


BaBar 2009 


25.45 ± 0.18 ± 0.28 


combined 


25.34 ± 0.06 ± 0.08 


combined 


25.20 ± 0.17 ± 0.28 



3. p — u> mixing 



In order to include the 1=0 contribution form uj — >• 7r + 7r~ we need to consider the corresponding symmetric 
(7, p, ui) 3x3 matrix propagator, with new entries Yl. yuj (q 2 ), H pcj (q 2 ) and q 2 — M 2 + H UJUJ (q 2 ), supplementing 
the inverse propagator matrix (23) by a 3rd row/column. Treating all off-diagonal elements as perturbations 
(after diagonalization) to linear order the new elements in the propagator read: 

n „ -IW) 



" (q 2 + n^(q 2 ))(q 2 -M 2 + n UJUJ (q 2 )) 



n ~ -iW<z 2 ) 



(q 2 - M 2 + n pp (q 2 )) (q 2 — M 2 + U^q 2 )) 



^"V-MS+IW^) ' (44) 

The self-energies again are the renormalized ones and in the two pion channel e + e~ — > tt + it~ given up to 
different coupling factors by the same self-energy functions as in the 7 — p sector. Thus, the bare self-energy 
functions read 

n 7 . = /V) , n p . = and n _ = fp f( q 2 ) , (45) 

and they are renormalized analogous to (14,15) subtracted at the oj mass shell. The p — uj mixing term 
is special here because if we diagonalize it on the p mass shell the matrix is no longer diagonal at the 
w-resonance, where 

n 2 a 2 — M 2 M 2 

n^V) - IW 9 2 ) - JL Rc U^(M 2 ) q -+ " U P UM 2 ) - -j-fl Re n pu (M 2 ) ± , (46) 
p p 

and which yields the leading 1=0 contribution to the pion form factor 1 . The uj induced terms contribute 
to the pion form factor 

^-^tt ^ (^) {Qlottit Quee) I^-)lo (.9 peedwinr ~l~ QueeQ pinr ) D pL0 \ J [e -D-y^J , (^^) 

which adds to (33). The direct e + e~ — > uj — > nir term given by — g^-^g^ee D uu by convention is taken 
into account as part of the complete w-resonance contribution. 

So far we have extended our effective Lagrangian by including direct p — uj, 7 — uj, ujtttt and wee vertices 
only, such that at the one-loop level only the previous pion loops show up. Missing are CL>7r + 7r~7r° and ojir 7 



"Typically, W™(M*) = Ml (h(M*) - Re h(Mf)) and D 7P (M,2) = - Sg? (^L^ggl . Similarly, 
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Fig. 9. Dynamical mixing parameter 5(E) obtained in our EFT, in contrast to the approximation by a constant. The latter 
seems justified by the narrow width of the u>. 



effective vertices, which are necessary in order to obtain the correct fuli w-width in place of the uj — > nir 
partial width only. Since the uj is very narrow we expect to obtain a good approximation if we use the proper 
full width in Im ILj^ = i M tJ T u (s), namely, 



r w -> r u (s) = r(w -> x, s ) = -^r w J Br ^ -> x ) 

X u { X 



Fx(s) 
Fx(Ml 



(48) 



where Br(V — > X) denotes the branching fraction for the channel X = 3tt,tt°^,2tt and Fx(s) is the phase 
space function for the corresponding channel normalized such that Fx (s) — > const for s — > oo [40] . 
If we include uj — p mixing in the usual way (see (1)) by writing 



e L> 77 + e (g p7r7r - g pee ) L> 7P - 5 pee g p7rTr D pp ■ ( 1 + 5j^BW u (s) 



/\e'D 7 J . (49) 



with BW w (s) — —M£/((s — M*) + iM w r w (s)) in our approach S e s(s) is given by 

c / \ peeQujiTTT QuieeQ ' pir-n') D puj 6 {^Quj^-n Qujee^) ^"fp /rn\ 

d ^y S ) = / TTS T)TTT / X I 50 ) 

which is well approximated by 

peeQujTTTr H~~ QuieeQ pinr ) ^-pu; (^) 



-Myn 



fj pTTTT fj pee 



Ml 



9 p-K-nQ pee 4o7T 

The second term g^eeQ 'pirn ~ 0.03 is an order of magnitude larger than than the first one (Jpee9wn-K ^ 0.003 
and thus is sensitive to g^ee once the gp^ has been fixed in the p-sector. In leading approximation 5 oc 
gujee/gpee ■ gpTr-n <?w^ir ■ The phase is actually fixed by the pion loop alone as we take couplings to be real 
(unitarity). We have \8\ = 1.945 x 10~ 3 and 4>& = 90.49°. 

A complete EFT treatment of the p — uj mixing, as well as the proper inclusion of the higher p's, requires 
the extension of our model, e.g. in the HLS version as performed in [13,14]. This is beyond the scope of the 
present study. Nevertheless, the discussion of the p — uj mixing presented above illustrates the need for a 
reconsideration of the subject. 
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4. Summary and Conclusions 

Our main point is to properly take into account the p — 7 mixing, which is responsible for the major part 
of the t vs. e + e~ discrepancy. The general message we have is that only a consequent application of the 
effective field theory approach can help to make progress in understanding low energy hadron data. Such 
attempts have been made recently within the HLS effective theory [13,14], where also a common consistent 
description (global HLS-model fit) of e+e~- and r-data has been found within a (p, to, (f>) mixing scheme. We 
have some difficulties to understand details of these elaborate calculations, but Fig. 8 of of [13] includes a 
component which has to be attributed to the p — 7 mixing, what we have been discussing here (our Fig. 6). 
Of course the effects considered here can easily be incorporated in other approaches like the ones proposed 
in [9] or [10]. 

We have based our "modeling" of the pion from factor on the low energy effective field theory of p, nn 
and 7, with the main assumption the pions to behave as point particles (sQED). We avoid some ad hoc 
elements, like imposing -FV(O) = 1 by hand, which is common practice when using GS like ansatze. Our 
result demonstrates two things: a) models should in any case be based on effective field theory, the "right" 
Lagrangian, to be sorted out by global fit strategies, b) obviously our simplest model has to be extended 
towards a full fledged resonance Lagrangians in order to be able to control the p — uj mixing and the energy 
range above the two kaon threshold. A proper treatment of the p — lo mixing requires an extension of the 
model to include the uj3it and ujttj couplings, which have anomalous parity, like 7r°77. Our study should 
be understood as a first step towards reconsidering the proper extraction of resonance parameters from 
experiments. 

In spite of the fact that we have to make use of a model which has its limitations, the relevant p — 7 
mixing effect (needed for example to correct the r data to be applicable for the evaluation of ajj ad ) can 
be determined from the e + e~-data solely. Hence, the r data represent independent additional information, 
which can be used to improve evaluations of hadronic effects which initially are directly related to e + e~-data. 

With this additional insight the original idea promoted in [41] indeed can work at the level of present 
standards in precision. Nevertheless, we should keep in mind that photon radiation effects from the composite 
hadrons are not fully under control and corresponding uncertainties are not easy to specify beyond the few 
per mil level. 

The model we use reproduces in some approximation the standard Gounaris-Sakurai model and improves 
it in several respects: as it should be we need no photon mass renormalization (which is intrinsically there in 
the GS model, as explained above), in addition to the s-dependence of the width (as incorporated properly 
in most versions of the GS model) we include the s-dependent p — 7 mixing, which leads to substantial 
modification of the GS model. The current conservation condition F v (0) = 1 is realized in our approach in 
a natural way (just by gauge invariance and by standard electromagnetic charge renormalization) as for the 
gauge invariant VMD type (7) and not by hand as it is done in the GS model, which is of the VMD type 
(6). 

What does it mean for the muon g — 2? It definitely shows that the r-decay isovector form factor must 
be corrected also for p — 7 mixing interferences, which means that relevant corrections not accounted for 
so far must be applied. These corrections are model dependent to some extent, as we assume pions to be 
point-like. Our calculation shows that the bulk of the effect is real and a 10% uncertainty seems to be a 
reasonable guess. The effects which only depend on the p parameters mass, width and leptonic branching 
fraction and for reasonable values of these parameters brings into fair agreement r-data based evaluations 
of ajj ad and the e + e _ -based ones. Thus phenomenologically, it reproduces rather precisely the pattern of 
the discrepancy between r and e + e~ extracted pion form factors (modulo differences which show up in the 
different measurements anyway). Our result strongly supports that the observed muon g — 2 discrepancy 
between theory and experiment is real and at the 3 a level (see e.g. [42], [38]). The r-data if properly corrected 
for isospin violating effects support this conclusion. 

For the lowest order hadronic vacuum polarization (VP) contribution to a M we find 

a had,LO[ e T ] = 690.96(1.06)(4.63) x lO" 10 (e + r) 
for the higher order vacuum polarization terms we find -206.68(0.36)(1.56) x 10" 11 , 103.89(0.16)(0.70) x 
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1CT 11 and 3.0(0.1) x 10" 11 , which adds to 

a had,VP,HO = _99. 79 (o.38)(0.86) x 10" 11 

(systematic errors of the first two essentially anti-correlated) . The corresponding updated value for the muon 
.9-2 is a^ c = 116591797(60) x lO" 11 which deviates from the experimental value ajf p = 116592080(54)(33) x 
10~ n [43] by aJfP - a^ llc = (283 ± 87) x lO" 11 corresponding to 3.3 a. 

It is important to note that the p — 7 mixing effect discussed here is included when evaluating the electro- 
magnetic current correlator, F 7T (s) and a^ ad from first principles in lattice QCD [44,45]. By comparison with 
the charged channel isovector correlator, the ratio F (s)/F^(s) could be measured within QCD, without 
reference to sQED or other hadronic models. 
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